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ABSTRACT 

The fundamental problem of how tunneling in thermal medium is completed is 
addressed, and a new time scale of order 1/friction for its termination, which is 
usually much shorter than the Hubble time, is pointed out. Enhanced non-linear 
resonance is responsible for this short time scale. This phenomenon occurs when the 
semiclassical periodic motion in a metastable potential well resonates with one of 
the environment harmonic oscillators coupled to its motion. 
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The usual picture of how the first order phase transition proceeds in cosmology, 
for instance the electroweak phase transition which is relevant for the scenario of 
electroweak baryogenesis is something like this; nucleation of the bubble of the 
true ground state (of broken symmetry) occurs with some probability in different 
parts of the universe. If the bubble formation does not take place strongly, the 
phase transition is terminated only when cosmological evolution changes the form of 
the potential (or more properly of the free energy) so that the symmetric phase is no 
longer the local minimum of the potential. This picture relies on the presumption 
that there is no inherent time scale for termination of the phase transition, or at least 
even if it exists, it is much larger than the Hubble time scale, which is the time for 
the potential change. We would like to exmaine this problem by working in a new 
formalism of the real-time description of the tunneling phenomena [§] , . Clearly, if 
the new time scale of the phase transition is much shorter than the Hubble time, one 
must reconsider the usual scenario. This might change the situation on, and even 
resurect the once failed GUT phase transition for the inflationary universe scenario 
i- 

Another important ingredient for our consideration is effect of cosmic environ- 
ment; the order parameter that describes the state of the universe such as the homo- 
geneous part of the Higgs field is necessarily coupled to matter fields that make up 
the thermal environment. One must then take into accout the presence of thermal 
medium in discussion of the tunneling rate and its time evolution, and estimate how 
dissipation due to the environment interaction modifies the basic tunneling rate. 

There are already many works on the subject of tunneling in medium 0, and 
most past works deal with a system in equilibrium as a whole. The Euclidean 
technique such the bounce solution is often used in this context 0, |^. Our 
approach here is different, and we attempt to clarify dynamics of time evolution 
starting from an arbitrary initial state of the tunneling system, which can be either 
a pure or a mixed state. We find it more illuminating to use a real-time formalism 
instead of the Euclidean method much employed in the literature. 

Our formalism is based on separation of a subsystem from thermal environment, 
and integrating out the environment variable including the interaction with the sub- 
system. In this picture dissipation seen in the behavior of the subsystem is due to 
our ignorance of huge environment degrees of freedom. Modeling of the environment 
and its interaction form with the subsystem is expected to be insensitive to the result 
one derives in this approach. We use the standard model of environment 0,0. 
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This method is best suited to a (by itself) non-equihbrium system which is im- 
mersed in a larger thermal equilibrium state. This methodology has been used by 



us in a number of problems in cosmology |T^, [|rT|. For instance, the relic abun- 
dance was calculated from this approach and the usual estimate of thermally 
averaged Boltzmann rate was justified at high temperatures. At very low temper- 
atures T <^ M (the mass of relic particle) our calculation differs from the usual 
one; the Boltzmann suppressed number density (MT/27r)^/^ e"^^/"^ at the freeze- 
out is replaced by some temperature power term. (There is some criticism against 
this calculation |12[.) But numerical importance of this effect is presumably minor, 
although it is a theoretically interesting issue. 

Here we consider the most basic problem of this kind, one dimensional subsystem 
described by a potential V{q). This system is put in thermal medium. The potential 
is assumed to have some local minimum at g = 0, which is separated at q = qb of the 
barrier top from a global minimum. The environment part is modeled by infinitely 
many, continuously distributed harmonic oscillators whose coordinates are Q{uj). Its 
coupling to the tunneling system is given by a Hamiltonian, 

q duj c{uj)Q{uj) . (1) 

The coupling strength is ciu) and Uc is some threshold frequency. Needless to say, 
one may imagine a generalized case in which the system variable q is the order 
parameter for the first order phase transition, the homogeneous Higgs field, and the 
environment oscillator Q{(jj) is a collection of various forms of matter fields coupled 
to the Higgs field. 

The basic equation in our problem is 

cPq dV f°° , / X d'^Q(uj) o ^, 

^ + — = -y^ duJc{Lo)Q{Lo), -^ + uj^Qiiu) = -ciiu)q. (2) 



One may eliminate the environment variable Q{uj) to get the Langevin equation [|13 



d'^q dV /■* , / x / x ^ / x 

where FQ{t) is linear in initial values of environment variables, Qi{uj) and Pi{uj); 

foo ( P(u;) \ 

FQ{t) = - / dujc{uj) Qi{uj)cos{ujt) + ^^sm{ut) . (4) 

By taking the thermal bath of temperature T = 1/(3 given by the density matrix, 

/ \ 1/2 

Pp{Q,Q') 



TT coth(/5cu/2) 



exp 



(g^ + g'^) cosh(/?cu) - 2gg' 



2 sinh(/?co' 

for each environment oscillator, the following correlation formula is obtained; 



(5) 



with 



({-^<3(t) 5-^g('5)}+)env = / (iu;r(u;) cosu;(r - s) coth(— ) = aR(r - s) , (6) 



r{uj) 



2uj 



(7) 



The kernel function aj in eq.@ is given by 

POO 

ai{t) = — / du r{u) sm{Lut) . (8) 

The combination, aR{t) + z ai{t), is a sum of the real-time thermal Green's function 
for g(ci;)— oscillators, added with the weight c^{uj). 

An often used simplification is the local, Ohmic approximation taking 
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with uJe = 0, which amounts to 



(9) 



(10) 
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a/(r) = 5iu^5{T)+ri5'{T). 
This gives the local version of Langevin equation, 

— + — + 6uj^ + — = 
dt"^ dq dt 

The parameter du"^ is interpreted as a potential renormalization or a mass renormal- 
ization in the field theory analogy, since by changing the bare frequency parameter 
to the renormalized uj\ the term 5uj'^ q is cancelled by a counter term in the poten- 
tial. On the other hand, rj is the Ohmic friction coefficient. This local approximation 



breaks down both at early and at late times [0 , but it is useful in many other cases. 
Any mixture of quantum states is described by a density matrix. 



(12) 



where |n) is a state vector for pure quantum state and < m„ < 1. The density 
matrix obeys the equation of motion; 



(13) 



where iJtot is the total Hamiltonina of the entire system. In configuration space this 
density matrix is given by its matrix elements, 



p(g,g'; , QV)) = , QMI P , Q'H) • (14) 

Its Fourier transform with respect to relative coordinates, q — q' , Q{uj) — Q'iyj) , is 
called the Wigner function and denoted by fw] 

/oo 
-oo 

ex^[-iip-i J duX{uj)P{uj)]p{q + i/2,q-i/2- Q{uj) + X{u) /2 ,Q{u) - X{uj) /2) . 

(15) 



It is easy to show from eq.(|T3D that this quantity obeys 

^fw dfw f J . I T-,t .\ ^fw . . 2 r\i \ '^fw 



p^- / duj P{uj)-^ + uj'Q{uj) 



dt dq J V ^^(^) 

(16) 



In one dimensional quantum mechanics the semiclassical approximation is excel- 
lent when the potential barrier is large, and we assume that this is also true in the 
presence of the system-environment interaction. In the semiclassical h ^ limit we 
have 

The resulting equation, being identical to the classical Liouville equation, has an 
obvious solution; 

fwiq,P,Q ,P) = J dqtdpi j dQidPi {qi ,Pi,Qi, Pi) 

■S {q - gel) S {p - Pel) S {Q - Qel) S {P - Pel) , (18) 

where qci{qi ,Pi ,Qi , Pi ;t) etc. are the solution of (||), taken as the classical equation 
with the specified initial condition. 

We consider the circumstance under which the tunneling system is initially in a 
state uncorrelated to the rest of environment. Thus we take the form of the density 
matrix, to get the reduced Wigner function after the trivial 
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Q{uj) , P{ijj) integration, 



fw\l^P'^t)= I dqidpif^' Jqi,Pi)K{q,p,qi,pi;t) 



(19) 



K{q,p, qi,Pi;t) = J dQidPi fw,Q{Qi , Pi) Hq - Qc\) 5 {p - pd) . (20) 

The problem of great interest is how further one can simplify the kernel function K 
here. 

In many situations one is interested in the tunneling probability when the envi- 
ronment temperature is low enough. At low temperatures of T ^ a typical frequency 
or curvature scale Ug of the potential V, one has 



o[Vt] < 



(21) 



In the electroweak and the GUT phase transition this corresponds to T ^ rriH {Higgs 
mass). Expansion of qd in terms of Qi{uj) , Pi{uj) is then justified. Thus, we use 



S{q-qci) = f — ^ exp [z Ag ( g - gc 
J Zn 



dXq 



exp 



(22) 



valid to the first order of Qi{uj) ,Pi{uj). A similar expansion for 5(p — Pci) using 
Pd\pci\uj),p[i\u), also holds. 

An alternative justification of the expansion (p2|) is to keep 0[c(co')] terms con- 
sistently in the exponent since both g^^'' (uj) and g^f ^ {u) are of this order. 

Gaussian integral for the variables Qi{uj) ,Pi{uj) first and then for Xg , Xp can 
be done explicitly with eg. (P^) , when one takes the thermal, hence Gaussian density 
matrix for the initial environment variable, pg'' . The result of this Gaussian integral 
leads to an integral transform of the Wigner function, f^^ (initial) (final), 
using the kernel function of 



Vd^ 

K[q ,p,qi,pi;t) = exp 

Ztt 



{q - qS\ P - Pd') J 



(0)^ 



(0) 

Q-Qci 

(0) 
P-Pcl 



(23) 



where the matrix elements of = are given by 



/ dui coth 


(5uj 


1 1 f 

— \Z[UJ 


t)\' 






UJ 


POD 

/ dto coth 




— \Z[UJ 




JlUc 




UJ 



'12 



'22 



1 r°° I3uj I 

duj coth — 9ft (z{uj , t)z*{uj , t)) 

2 UJ 



hi 
21 



(24) 
(25) 
(26) 
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where z{uj ,t) = q''^\uj ,t) +iuj q''J'\uJ ,t) , and z{uj ,t) = p^^\uj ,t) + iu pli\uj ,t) . 

Quantities that appear in the integral transform are determined by solving differ- 
ential equations; the homogeneous Langevin equation for q''^^ and an inhomogeneous 
linear equation for z{u ,t) and z{u ,t), 

^(dy\ ^ 2 rf. aj{t - s)qf{s) = , (27) 



dt'^ \dq J (0) ^0 



cl 



d^4^ ,t) dHr\ ^^^^^ rt ^^^^ _ ^^^^^ ^ s) = - c{uj)e'^' . {2i 



df^ \ dq^ I (0) ' JO 



cl 



A similar equation as for z{u! ,t) holds for z{uj ,t) . The initial condition is 
q^^\t = 0) = q,, p(J)(t = 0) =p,, z{u;,t = 0) = 0, z{u; ,t = 0) = . 

The physical picture underlying the formula for the integral transform, eq.(|l9|) 
along with (p3D , should be evident; the probability at a phase space point {q,p) is 
dominated by the semiclassical trajectory q^^^ (environment effect of dissipation being 
included in its determination by eq.(p7D) reaching {q ,p) from an initial point (g^ ,pi) 
whose contribution is weighed by the quantum mechanical probability f^) initially 
given. The contributing trajectory is broadened by the environment interaction 
with the width factor Jlij- The quantity In given by (p^, for instance, is equal to 



(g — g^i*^)^ ^ environment driven fluctuation under the stochastic force FQ{t). 

We note that in the exactly solvable model of inverted harmonic oscillator po- 
tential the identical form of the integral transform was derived P] without resort to 
the semiclassical approximation. Explicit form of the classical and the fluctuation 
functions {q^^^ and z{uj ,t)) is given there. 

To proceed further, one may separate the tunneling potential V{q) at the barrier 
top location, q = qs- We distinguish two cases of the potential type, depending on 
the value of V{oo) relative to the local minimum value Vm at g = in the potential 
well. When ^(oo) < Vm, the classical motion in the overbarrier region g > is 
monotonic ending at g = oo, while for V{oo) > Vm the motion is a damped oscillation 
towards go at the true minimum, unless the friction is large. If the friction is larger 
than a critical value of ~ 2uj^, with cj* being the curvature of the potential at the 
global minimum, there occurs the overdamping such that g^°^ go monotonically. A 
typical interesting case for V{oo) > Vm is the asymmetric double well as may occur 
in the first order electroweak phase transition 

We consider here a few problems to illustrate consequences of our general formula 
of the integral transform. One problem is calculation of the barrier penetration fac- 
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tor; to determine the outgoing flux in the overbarrier region for the type of potential 
of V{oc) < Vm, assuming an initial energy eigenstate under the potential V. The 
other problem is the tunneling probability for the asymmetric wine bottle type of 
potential. 

Energy eigenstates are special among pure quantum states, since they evolve only 
with phase factors. Thus, if one takes for in eq.(|T2|) the energy eigenstate of 

the total Hamiltonian, the density matrix does not change with time. On the other 
hand, if one takes the energy eigenstate for the subsystem Hamiltonian, then the 
density matrix changes solely due to the environment interaction. It is thus best 
to use a pure eigenstate, or its superposition, of the subsystem Hamiltonian when 
one wishes to determine how the barrier penetration factor is modified in thermal 
medium. 

We first discuss the stationary phase approximation for the initial state. Starting 
from the formula. 



J di exp 



-ipii + \Yiil)l{qi - |) + lnV'„(gi + |; 



(29) 



we locate the stationary point by ^ = 0,^ = for each exponent factor. The par- 
tial derivative ^ = is taken with the understanding that the rest of pi integration 
contains smooth functions of pi. This leads to the stationary point at 

e = 0, p^ = In{q^) (30) 

"^"^^ " ^ WM' ■ ^ ^ 

Here /„(x) is the usual flux factor at x for the pure state \n). This gives the result 
for the initial density matrix, 

f^^{q^,P^)-Y.^n\i^n{qi)?2'K5{p,-IM) . (32) 

n 

The barrier penetration factor is calculated from the flux formula, 

i{q.t)= r ^pfi?\q,p;t). (33) 

Thus, the semiclassical plus the stationary phase approximation gives 

iq-qS^r' 



I{q,t)^ dqi\ilj{q^)\ W — — exp 

Jq, V ^^-'11 



21 
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(34) 
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for a pure initial state given by the wave function ip{x), where is the turning point 
in the overbarrier region. 

In the region of g ^ there is always a classical trajectory that reaches the 
point q of the Gaussian peak from an initial qi in the range qi > q^, and the entire 
region within the Gaussian width \/Tii is fully covered by the integral. The factor 
outside the Gaussian exponent is well approximated by its peak value in the weak 
coupling case. Thus, one obtains 

J(g,t) ^ |^(a;o)r (-i^o) = mx,)\' {x, {^-^^^^ , (35) 

where Xo{q,t) is determined from 

qS\'li = XQ,Pi = I{xo),Qiiuj) = 0,Pi{uj) = 0;t) =q. (36) 
We now take the WKB wave function in the overbarrier region, 



T{E) r ri^ 



i^iqi) = r-— exp 
'P{(li) 



dxp{x) 



p{x) = ^2{E-V{x)) , (37) 



to get a factorized form of the flux 



f^E^[^) . ,38, 

A feature that characterizes the classical trajectory q^^^ is of special interest; its 
initial energy is 

H,{t = 0)^^piq,f + Viq,) = E. (39) 

This formula for the flux reproduces our previous result for a specific potential 
of the inverted harmonic oscillator 0, V{x) + ^Su^x"^ = —■^uj^x'^ , with uji the 
renormalized curvature. The present approach actually improves our previous result; 

^ ^ 9 xo + gljxo) _^ 9 +ujBg 
gI{xo) + guj^xo usig + ujb9) ' 



where ujb ~ u;^ — 77/2 is a diagonalized frequency. The limiting formula of eq.( 
with /(a;o) ojbXq as xq 00, is valid in the infinite q limit, as derived in ||^. The 
function g{t) is the homogeneous solution of the Langevin equation given explicitly 
ini; 

J\[2 poo 

g{t) = — sinh{ujBt)+2 / du H (u) sm{ujt) (41) 

LOB Ji^c 

^('^) ^ , 2^ 2^21^ , ,,2 ' (42) 

{u^ + uj%y + {■7Tr{uj)y 

/•OO 

iV^ = 1 - 2 / duH{uo)uj. (43) 



Both at early and late times the factor / ~ 1 , deviating from unity only for the time 
range of order I/ojb- 

According to the view of the potential V{q) that determines the semiclassical 
penetration factor |T(i?)p should be expressed in terms of the renormalized parame- 
ters, ojR in this case. This explains the bulk of the suppression caused by dissipative 
environment interaction, as explained in [Q. 

For discussion of a more general case of finite V{oo) < Kn we use the local, Ohmic 
approximation, which becomes excellent at late times. A potential that decreases 
fast as g — s> oo is assumed; ^ — > . The acceleration term can then be neglected 
when the friction satisfies rf ^ This is a slow rolling approximation, and it 

always holds for q large enough. The classical equation is then solved as 

,dV.. 



V 



/■<? dV 1 



which gives the factor 



f-^-— ^ (45) 
dq^ r]p{oo) 

Thus, the tunneling probability decreases with time along with the decreasing slope 
of the potential. This result poses a curious question; the tunneling may not be 
terminated within a limited finite time. We shall encounter a similar situation for 
the asymmetric wine bottle potential. 

We next consider the case in which the potential is very steep at both ends; 
V^(±oo) = oo . The tunneling rate, from the inner region at q < qs into the outer 
region at g > g^, is an important measure of tunneling phenomena and is given by 
the flux at g = g^; P{t) = — I{qB , t) , which is equal to 



(qb - ge l 

2/ii 



(0)^2 



(46) 



On the other hand, the tunneling probability into the overbarrier region at g > x is 
given by 



/(■\ 1 
dqidpi f^' {qi ,pi) du ^— ^ exp 



{u - q^^ f 



21 
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• (47) 



In both of the quantities, P(t) and P{qB ,t) (the tunneling probability into g > qs) 
it is essential to estimate how q^^^ and In varies with time. 
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We consider an initial state localized in the potential well so that the dominant 
contribution in the {qi,Pi) phase space integration is restricted to qi < qb- In the 
rest of discussion anharmonic terms play important roles, but at first we work out 
the harmonic approximation; 

V{q) + ^Su;'q'^^u;U\ (48) 

near the bottom of the well at q = 0. In the Ohmic approximation the classical 
solution and the fluctuation is given by 



q^^i' = cos uot + sin Co^t e"''*/^ q, + — ^ e"''*/^ , (49) 

, X c(a;) 
z{u:,t) = - V • 

UJ^ — UJq— ILOT] 

V 2cuo 2cuo / ' 

using 0)0 = \J^o ~ ^ ■ III f rest of discussion we assume a small friction, 77 ^ a;o . 
Near uj — ujq the fluctuation is approximately 

Z{U , t) « ( ^ ^iC.,t _ 1 ~ ^ A g-W2 . (51) 

^ ' ^ a; + cuo-^r7/2 V (Dq / ^ ^ 

This formula is valid at t < Ofl/r^]. The appearance of the linear t term is a resonance 
effect. The resonance roughly contributes to h\{t) by the amount, ryi^e"''* x a 
smooth uj integral which is cut off by a physical frequency scale. Thus, the width 
factor ill initially increases with time until the time scale of order I/77. 
The width factor asymptotically behaves as 

/ii(i) = /ii(oo) + 0[e-''*/2], (52) 
7^^(00) f« — + — + -\t\ (53) 

The asymptotic value of /ii(oo) has the familiar zero point fluctuation of harmonic 
oscillator and in the last term the dominant finite temperature correction, valid for 
this Ohmic model at T <^ a;o. In any event the probabihty rate P{t) finally decreases 
to zero, along with 

pf + ^{lB-,f)^0. (64) 

Mil 

Moreover, the final tunneling probability P{qB , 00) has a finite value, and typically 
is very small for a large potential barrier. For instance, for the asymmetric wine 
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bottle potential shortly discussed, 



P(gB,oo)^i /-^e-«^^>o, (55) 
4 y ZnVh 

with Vh the barrier height much smaller than c^o- This poses again a curious question; 
it appears that decay of a prepared metastable state localized in the potential well 
is never completed. 

This simple picture is however valid only when one ignores anharmonic terms in 
the tunneling potential, but they must be there in order to give any realistic tunneling 
potential. The most important in the following discussion is effect of anharmonic 
terms in the equation for the fluctuation z{uj ,t). Presence of anharmonic terms 
gives a non-trivial periodicity in the coefficient function (^r) (o) for (PS]), assuming 

^cl 

a small friction r] <^ cuq. 

The homogeneous part of the z{uj ,t) solution which might exhibit the well known 
parametric resonance is closely related to the behavior for our inhomogeneous 
z{uj , t) solution. It is then important to check whether the relevant parameter in the 
periodic coefficient function (o) falls in the instability or the stability band. As 

^cl 

is shown in [0, unbounded exponential growth of does not take place. On the 
other hand, the power-law growth is observed in numerical computation. Moreover, 
we find that our z{uj ,t) solution belongs to the boundary between stability and 
instability bands. At the resonance frequency the enhancement factor due to the 
boundary effect is much larger than what one might expect from the harmonic case 



(|5lD. We shall interpret this phenomenon as influenced in a subtle way by the 
parametric resonance, although it is not the parametric resonance itself. 

It is best to discuss the resonance enhanced tunneling mechanism in concrete 
examples. We take the asymmetric wine bottle potential, which is described in the 
well and its vicinity region by 

V{q)-^{q'-2qBqf. (56) 

The curvature parameters at two extema of g = and q = qs are Uq = 2Ag^ , cu^ = 
Ag^ , and the barrier height seen from the bottom of the well is \4 = = uj^q^/S . 
The classical q^^\ and the fluctuation z{ul) ,t) equations in the Ohmic approximation 
are written using rescaled variables, y = q^i^ /qs and r = uJot/2, 

y" + y(y-l)(y-2) + ^y' = 0, , (57) 



z"+{4-6{y--y'))z + ^z' = - -A^ ^../.o ^ (5^ 
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where ' = d/dr. 

For a small friction one has approximate forms of solution in terms of the Jacobi's 
elliptic function; 



(y - 1)2 = 1 + - 2v^ sn^ 1^ Vl + v^ ( ^ + r 



du 



p^°^(t) = 2i/2^sn vTTVe ( ^ + const ) , A; j cn v^l + ( ^ + const ) , fc^ 



1 



2^t 



(59) 



(60) 
(61) 



These formulas are valid for t <^ l/r] . There are resonances at uq , 2ljo , ■ ■ ■ nuo , ■ • •, 
where 




du 



1-m2)(1- 



(62) 



We found so far that these explicit solutions are not very illuminating, and nu- 
merically integrated the coupled q'^^ and z{uj ,t) equations. In the u integral (|^) for 



the width factor In the largest contribution is found to come from the fundamen- 
tal resonance at u = ujq, then contribution from higher harmonics at nCoo follows. 
Thus, phenomenon of a non-linear resonance occurs. At a time of Ofl/r/] the width 
factor III becomes maximal at its value much larger than in the harmonic case. The 
decrease observed at late times is due to the friction; we explicitly checked that for 
the zero friction, the fluctuation \z{uj ,t)/c{uj)\'^ increases in time without a bound, 
with an averaged time power close to 4 at the resonance. Note that even if the effect 
of the friction is turned off, there exists an important environment effect here; the 
environment interaction drives the non-linear resonance oscillation. We refer to our 
paper [|16] for detailed numerical results. 

We numerically checked [|I^ the behavior of the most important part given by 
the product of two competing exponential factors. 



A 



exp 



— tanh 



(3uJo + ujI qf 



X exp 



JO) )2 



21 
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(63) 



the one for the initial state and the other for the kernel factor in the probability 
rate. Remarkably, the largest contribution comes, not from the dominant initial 
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component near the zero point energy, rather from the initially suppressed excited 



component. The first exponent in (|63|) is in proportion to Ei, while the second one 
goes roughly like E~'^, hence a maximum may appear somewhere away from the 
lowest energy of Ei = ujq/2. In the example of t]/ujq = 0.0025 the maximum product 
factor is of order 10~^ at Ei ^ Ax ujo/^, 6 orders of magnitudes larger than what one 
expects from the lowest energy state and also the asymptotic value of order 10~^^. 
A large value of the product factor of order unity suggests an interesting possibility 
of a rapid and violent termination of the tunneling. 

The reason why one observes enhanced tunneling at resonant frequencies of 
UJ = nuQ is as follows. The semiclassical g— motion in the left potential well has 
a periodicity 2tt/Cjq if one neglects the friction. The environment oscillators act 
as stochastic fiuctuation to this motion in such a way that the effective potential 
including one environment oscillator, 

V{q) + qc{uj) (^Qi{u) cos{ujt) + '^^^ sin{ujt)^ , (64) 

changes the potential barrier periodically. When the g— particle hits against the 
potential wall, the potential barrier may become smallest, hence the tunneling rate 
maximally enhanced, if one of environment oscillators has a period exactly equal 
to that of the classical motion, uj = huq. This is precisely the condition of the 
non-linear resonance we have been discussing. The resonance enhanced tunneling 



thus envisaged seems to have little connection to the stochastic resonance much 
discussed in the literature. 



In summary, we gave a new mechanism of resonance enhanced tunneling, using 
the real-time semiclassical formalism. Our approach suggests a new time scale of 
order l/t], the inverse friction, for completion of the first order phase transition. 
Evidently much has to be done to apply the idea here to realistic problems. 

I wish to thank my collaborator, Sh. Matsumoto who shared most of the results 
presented here. 
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